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Abstract 

Heavy resonances appearing in the clean Drell-Yan channel may be the first new physics 
to be observed at the proton-proton CERN LHC. If a new resonance is discovered at the 
LHC as a peak in the dilepton invariant mass distribution, the characterization of its spin 
and couplings will proceed via measuring production rates and angular distributions of the 
decay products. We discuss the discrimination of the spin-1 of Z' representative models 
(Zggp^^, Z^, Z^, Z^, Zlpj^, and Z^lr) against the Randall- Sundrum graviton resonance 
(spin-2) and a spin-0 resonance (sneutrino) with the same mass and producing the same 
number of events under the observed peak. To assess the range of the Z' mass where 
the spin determination can be performed to a given confidence level, we focus on the 
angular distributions of the Drell-Yan leptons, in particular we use as a basic observable 
an angular-integrated center-edge asymmetry, Ace- The spin of a heavy Z' gauge boson 
can be established with Ace up to Mz' — 3.0 TeV, for an integrated luminosity of 100 
fb~^, or minimal number of events around 110. We also examine the distinguishability of 
the considered Z' models from one another, once the spin-1 has been established, using 
the total dilepton production cross section. With some assumption, one might be able to 
distinguish among these Z' models at 95% C. L. up to Mz' ~ 2.1 TeV. 
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1 Introduction 



Heavy resonances with mass around 1 TeV or higher are predicted by numerous New 
Physics (NP) scenarios, candidate solutions of conceptual problems of the standard model 
(SM). In particular, this is the case of models of gravity with extra spatial dimensions, 
grand-unified theories (GUT), electroweak models with extended spontaneously broken 
gauge symmetry, and supersymmetric (SUSY) theories with i?-parity breaking (-^p). 
These new heavy objects, or 'resonances', with mass M 3> M\y^Zi may be either pro- 
duced or exchanged in reactions among SM particles at high energy colliders such as the 
LHC and the International electron-positron linear collider (ILC). A particularly inter- 
esting process to be studied in this regard at the LHC is the Drell-Yan (DY) dilepton 
production (/ = e, /i) 

p + p^l^l- +X, (1.1) 

where exchanges of the new particles can occur and manifect themselves as peaks in 
the invariant mass M. Once the heavy resonance is discovered at some M = Mr, 

further analysis is needed to identify the theoretical framework for NP to which it belongs. 
Correspondingly, for any NP model, one defines as identification reach the upper limit 
for the resonance mass range where it can be identified as the source of the resonance, 
against the other, potentially competitor scenarios, that can give a peak with the same 
mass and same number of events under the peak. This should be compared to the discovery 
reach, which specifies the (naturally more extended) mass range where the peak in the 
cross section pertaining to the model can just be observed experimentally. Clearly, the 
determination of the spin of the resonance represents an important aspect of the selection 
among different classes of non-standard interactions giving rise to the observed peak. 

Tests of the spin-2 of the Randall-Sundrum [1] graviton excitation (RS) exchange in 
the process (11.11) at LHC, against the spin-1 hypothesis, have been recently performed, 
e.g., in Refs. [2-4] on the basis of the lepton differential polar angle distribution, and in 
Ref. [5] using the azimuthal angular dependence. In the reverse, the identification of the 
spin-1 Z's has been discussed in [6,7]. The above-mentioned differential angular analysis 
in the polar angle has been applied to the search for spin-2, spin-1 and spin-0 exchanges 
in the experimental studies of process (11.11) at the Fermilab Tevatron proton-antiproton 
collider [8]. 

In Ref. [9], the discrimination reach at the LHC on the spin-2 RS graviton resonance 
or, more precisely, the simultaneous rejection of both the spin-1 and spin-0 hypotheses for 
the peak, has been assessed by using as basic observable an angular-integrated center-edge 
asymmetry, Ace, instead of the 'absolute' lepton differential angular distribution. The 
potential advantages of the asymmetry Ace to discriminate the spin-2 graviton resonance 
against the spin-1 hypothesis were discussed in Refs. [10,11]. 

Here, along the lines of Ref. [9] but in the reverse direction, we apply the same basic 
observable Ace, to the spin-1 identification of a peak observed in the dilepton mass distri- 
bution of process (II. ip at the LHC, against the spin-2 and spin-0 alternative hypotheses. 
For explicit NP realizations, for the spin-1 Z' models we refer to Refs. [12]; for the alter- 
native spin-2 and spin-0 hypotheses we refer for the RS graviton resonance to [1] and for 
the SUSY sneutrino exchange to [13,14], respectively. 

It turns out that Ace should provide a robust spin diagnostic for the spin-1 case 
also. Moreover, we examine the possibility, once the spin-1 for the discovered peak is 
established, of differentiating the various representative Z' models from one another. For 
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this purpose, we must use the total dilepton production cross section or, equivalently, the 
rate of events of reaction fll.ll] under the peak. Identification of Z' models have been 
discussed recently in, e.g. [7, 15, 16] with different sets of observables, namely, forward- 
backward asymmetry Afb on and off the Z' resonance, Z' rapidity distribution, cross 
section times total width, a x F^/, as well as in different processes [17,18]. It was found 
that, on the basis of Afb only, pairs of Z' models become indistinguishable at a given level 
of significance, starting from relatively low values of M^/ of the order of 1-2 TeV, even 
at £int much higher than 100 fb~^. These ambiguities can be reduced by the combined 
analysis of the observables mentioned above, and at = 100 fb~ , some models could 
be discriminated up to Z' mass of the order of 2-2.5 TeV. As we will note below, on the 
basis of a simple criterion, the precise determination of the total cross section itself 
might provide a somewhat stronger discrimination potential, in the sense that all models 
could be pairwise distinguished from one another up to Z' masses of about 2.1 TeV. 

In Sec. 12] we present a brief introduction to the main features of the different models 
considered in the analysis, and the expected relevant statistics; Sec. |3] is devoted to the 
spin-1 identification of Z' bosons against the spin-2 RS and the spin-0 sneutrino hypothe- 
ses; in Sec. S] we derive the differentiation of Z' models among themselves obtainable 
at the LHC from consideration of total dilepton cross sections, wheras Sec. E] is devoted 
to a brief discussion of the reduced-energy, low-luminosity domain relevant to the early 
running period of the collider. Finally, Sec. [6] contains some conclusive remarks. 



2 Cross sections and considered NP models 

For completeness and to fix the notations, we start by recalling the basic expression for 
the cross section of process (11. ip . and present a mini-review of the NP models we want 
to compare. 

The parton model cross section for inclusive production of a dilepton with invariant 
mass M can be written as 

da{Ru) -j,^V/.(6,M)/,(6,M)^(.+,^ /- + /-). (2.1) 



Here, s is the proton-proton center-of-mass energy squared; z = cos^c.m. with ^c.m. 
the lepton-quark angle in the dilepton center-of-mass frame; y is the dilepton rapid- 
ity; fij{^i^2, M) are parton distribution functions in the protons Pi and P2, respectively, 
with ^1^2 = {M/ s/s) exp(±|/) the parton fractional momenta; finally, ddij are the partonic 
differential cross sections. In (12.11) . the factor K accounts for next-to- leading order QCD 
contributions [19,20]. For simplicity, and to make our procedure more transparent, we 
will use as an approximation a global fiat value K = 1.3. 

Since we are interested in a (narrow) peak production and subsequent decay into 
the DY pair, pp R ^ we consider the lepton differential angular distribution, 

integrated over an interval of M around Mr: 

da{Ru) ^^r^ da ^ 



Jmr~aai/2 J^Y^Mdydz 
The number of events under the peak, that determines the statistics, is therefore given 
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by: 



a{Ru) = a{pp ^ R) ■ BR{R ^ /+/-) = dz dM dy———. (2.3) 

J-z.ut J Mr- AM/2 J-Y dMdydz 

For the full final phase space, Zcut = 1 and Y = \og{^/s/M). However, if the finite 
detector angular acceptance is accounted for, Zcut < 1 and Y in Eqs. (12.21) and (12.31) must 
be replaced by a maximum value ymax{z, M). Concerning the size of the bin AM, it should 
include a number (at least one) of peak widths to enhance the probability to pick up the 
resonance. In the models we will consider, widths are predicted to be small, typically of 
the order of a percent (or less) of the mass Mr, so that the integral under the peak should 
practically be insensitive to the actual value of AM. Conversely, the SM 'background' is 
expected to depend on AM. In our analysis, we adopt the parametrization of AM vs. 
M proposed in Ref. [6] and, denoting by Nb and Ns the number of 'background' and 
'signal' events in the bin, the criterion Ns = 5\/Nb or 10 events, whichever is larger, as 
the minimum signal for the peak discovery. 

To evaluate the statistics, we shall use in Eqs. (12. 2p and (12.31) the CTEQ6.5 parton 
distributions [21], and impose cuts relevant to the LHC detectors, namely: pseudorapidity 
|?7| < 2.5 for both leptons assumed massless (this leads to a boost-dependent cut on z [11]); 
lepton transverse momentum p±_ > 20 GeV. Moreover, the reconstruction efficiency is 
taken to be 90% for both electrons and muons [22] and throughout this paper, except for 
Sec. [5l a time-integrated LHC luminosity Cmt = 100 fb~^. 

For the proton-proton initiated process (II. ip . only the 2;-even parts of the partonic 
differential cross sections contribute to the right-side of Eq. (12.21) . 2;-odd terms do not 
contribute after the |/-integrationj3 Also, due to Mz <C M/j and the narrow width peak, 
the resonant amplitude interference with the SM is expected to give negligible contribu- 
tions to the right-hand sides of (12. 2p and (12. 3p after the symmetric M-integration around 
Mr needed there. Thus, we can retain in these equations just the SM and the resonance 
pole contributions!^ This fact was noticed for Z'-exchange in, e.g., Refs. [6,19], but holds 
also for the RS graviton and for the scalar sneutrino exchanges discussed later. 

2.1 Z' models 

In a wide variety of electroweak theories, in particular those based on extended, spon- 
taneously broken, gauge symmetries, the existence of one (or more) new neutral gauge 
bosons Z' is envisaged. These additional gauge bosons could be accessible at the LHC. 
A new neutral gauge boson would induce additional neutral current interactions. The 
color-averaged differential cross section for the relevant, leading order, partonic subpro- 
cess qq ^ Z' ^ can be expressed asjfl 

^Accordingly, for the qq and gg subprocesses, only the combinations of parton distributions 
[/,(Ci, M)/,-(6, M) + fq{ii,M)fq{^2, M)] and M)/s(6, M) arc effective in the cross sections 

and ([Q]) . 

^Actually, such interference can in principle contribute appreciably to the differential cross section 
da/dAfdy [9], and plays a role in the forward-backward asymmetry (which we do not consider here). 
•^We neglect fermion masses as well as potential effects from the (tiny) Z — Z' mixing. 
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with 



•Z' 



(2.5) 



Xz' = 



Eq. fl2.4l) shows that the spin-1, Z', exchange in process f 1 1.11) has the same symmetrical 
angular dependence as the SM 7 and Z exchanges. 

The list of Z' models that will be considered in our analysis is the following: 

(i) The three possible U{\) Z' scenarios originating from the exceptional group 
spontaneous breaking. They are defined in terms of a mixing angle /3, and the cou- 
plings are as in Tab. [H The specific values = 0, = 7r/2 and /3 = arctan — -\/5/3, 
correspond to different breaking patterns and define the popular scenarios Z^, 
Z'^ and Z'^, respectively. 

(ii) The left-right models, originating from the breaking of an 5*0(10) grand- unification 
symmetry, and where the corresponding Z[^yi couples to a combination of right- 
handed and B — L neutral currents [B and L denote lepton and baryon currents), 
specified by a real parameter olr bounded by a/2/3 ^ olr ^ Corresponding 
Z' couplings are reported in Tab. [H We fix olr = v^, which corresponds to a pure 
L-R symmetric model. 

(iii) The Z^lj^ predicted by the 'alternative' left-right scenario. 

(iv) The so-called sequential Zggj^, where the couplings to fermions are the same as those 
of the SM Z. 

Detailed descriptions of these models, as well as the specific references, can be found, e. g., 
in Ref. [12]. All numerical values of the Z' couplings needed in Eq. (12. 5p are collected in 
Table [H where: A = cos (3/2^ and B = v^sin/3/12 are used. 

Current Z' mass limits, from the Fermilab Tevatron collider, are in the range 800 — 
1000 GeV, depending on the model [23]. 

2.2 RS graviton excitation 

We consider the simplest scenario in the class of models based on one compactified warped 
extra dimension and two branes, proposed in the context of the SM gauge-hierarchy 
problem in [1]. The model predicts a tower of narrow Kaluza-Klein (KK), spin-2, graviton 
excitations G^"-* (n > 1) with the peculiar mass spectrum M^*^^ = M^^^Xn/xi {xi are the 
zeros of the Bessel function, Ji{xi) = 0). Their masses and couplings to the SM particles 
are proportional to A^^ and 1/A^, respectively, with A^ the gravity effective mass scale on 
the SM brane. For A^ of the TeV order, such RS graviton resonances can be exchanged 
in the process (11. ip and mimic Z' exchange. The independent parameters of the model 
can be chosen as the dimensionless ratio c = k/Mp\ (with k the 5-dimensional curvature 
and Mpi = I/v^SttGn the reduced Planck mass), and the mass Mq of the lowest KK 
resonance G^-^\ Accordingly, A^ = Mc/cxi. 

The differential cross sections for the relevant partonic subprocesses needed in (12. 2p 
and (O, qq^G ^ l+l' and gg ^ G ^ /+/", read, with k = v^cxi/M^ [2,24-27] 
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[Agg{z) + Ag,{z)] \xg 
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Xg = 



M2 
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Table 1: Left- and right-handed couphngs of the first generation of SM fermions to the 
Z' gauge bosons, in units of l/c^y for the and LR models, and l/svi/Cvi/\/l ~ 2s^ for 
the ALR model [9], where cw = cos^iy, sw = sin6'iy. 
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where 

The theoretically 'natural' ranges for the RS model parameters are 0.01 < c < 0.1 and 
A^ < lOTeV [28]. Current lower bounds at 95% C.L. from the Fermilab Tevatron collider 
are, for the first graviton mass: Mq > 300 GeV for c = 0.01 and Mq > 900 GeV for 
c = 0.1 [23,29]. 



2.3 Sneutrino exchange 

Sneutrino (t>) exchange can occur in SUSY with i?-parity breaking, and represents a 
possible, spin-0, interpretation of a peak in the dilepton invariant mass distribution of the 
process (11.11) . The cross section for the relevant partonic process, — >• z/ — > is flat 
in z and reads [13]: 



dz N,AM^\e^ \Xu\ oqd, X. M^-Ml + iM^V^ 



In Eq. (12. 8p . A and A' are the i?-parity- violating sneutrino couplings to and dd, 

respectively. Actually, in the narrow-width approximation, the cross section (12.81) turns 
out to depend on the product X = [X'^Bi, with Bi the sneutrino leptonic branching ratio. 
Current limits on X are rather loose [30], and we may consider for this parameter the 
range 10"^ <X< 10"^. For 10"^ <X< 10"^ the range is > 280 - 800 GeV 



2.4 Statistics and model signature spaces 

In Fig. [1], we show the predicted number of resonance (signal) events Ns in the Drell- 
Yan process (II. ip at LHC, vs. Mr, where R = Z', G, v denotes the three alternative 



6 




1 2 3 4 5 

Mr (TeV) 

Figure 1: Expected number of resonance (signal) events Ns vs. (i? = Z', G, 9^) at the 
LHC with £i„t = 100 fb"^ for the process pp ^ R ^ l^l + X {I = e,fi). Event rates 
for various Z' models are shown. Green area corresponds to graviton signature space for 
0.01 < c < 0.1 while the yellow area is the sneutrino signature space for 10~^ < X < 
10"^. Minimum number of signal events needed to detect the resonance (S-cx level) above 
the background and the minimum number of events to exclude the spin-2 and spin-0 
hypotheses at 95% CL. are shown. Error bars correspond to the statistical uncertainties 
for the ALR model. 

possibilities outlined in the previous subsections. The assumed integrated luminosity is 
Cint = 100 fb~^, the cuts in phase space relevant to the foreseen detector acceptance 
specified above have been imposed, and the channels / = e, /i have been combined. Also, 
the minimum signal for resonance discovery above the 'background' at 5a is represented 
by the long-dashed line. 

For any model, one can define a corresponding signature space as the region, in the 
(Mr, Ns) plot of Fig. [H that can be 'populated' by the model by varying its parameters 
in the domains mentioned above. Clearly, in regions where the signature spaces overlap, 
the values of Mr are such that it is not possible to distinguish a model as the source of 
the peak against the others, because the number of signal events under the peak can be 
the same. Further analyses are needed in these cases to perform the identification of the 
peak source. For example, the 'blue' area in Fig. [1] corresponds to the graviton signature 
space for 0.01 < c < 0.1, while the yellow area (which has substantial overlap with the 
blue one — indicated as green) is that for the sneutrino signature space corresponding to 
10-5 < X < io-\ 
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As regards the discovery and identification of Z' we are interested in, tlie signature 
spaces in Fig. [T] reduce to tlie lines labelled by the different models, because the event 
rates are fixed, once is given, through the couplings in Table 1. Fig. [T] shows that, 
with the assumed (high) luminosity of 100 fb~ , Z' gauge boson masses up to 4-5 TeV 
are in principle within the 5-0" reach of the LHC, consistent with earlier studies [12]. We 
here assume that the Z' can only decay to pairs of SM fermions in order to obtain the 
leptonic branching ratio Bi. It is important to note that in many models, where Z' can 
also decay to exotic fermions and/or SUSY particles this overestimates Bi and, thus, the 
search reach. 

On the other hand. Fig. [1] demonstrates that, as far as the production rate of DY pairs 
is concerned, there is a substantial overlap between the Z' and the v signature spaces, 
which determines a domain in (Mj>, X) where spin-0 v exchange and Z' exchanges are not 
distinguishable because they lead to the same event rate under the peak. The same is 
true for the spin-2, RS model. However, as shown by Fig. [H in this case it is interesting 
to note that, if one literally takes the suggested range c < 0.1 as the 'naturally' preferred 
one, the ALR and SSM scenarios can be discriminated against the RS (spin-2) resonance 
already at the level of event rates in a wide range of mass values accessible to the LHC, 
with no need for further analyses. Conversely, only the and LR Z' models possess a 
'confusion region' with the RS resonance G, concentrated near the upper border of the 
graviton signature domain. 



3 Identification of Z' spin-1 

We now turn to the identification of the spin-1 of the Z' boson vs. the spin-0 and spin-2 
hypotheses using the angular distribution of the final-state leptons. 

For this purpose, we adopt the integrated center-edge asymmetry Ace, defined as 
[10,11]: 



ct{Rii) 



with acE{Ru 



dajRu 
dz 



■dz. (3.1) 



In Eq. f l3.ll) . R denotes the three hypotheses for the resonances we want to compare, 
namely: spin-1 (V); spin-2 (G); and spin-0 (S). Moreover, < z* < Zcut is an a priori 
free value of cos6'cm that defines the 'center' and 'edge' angular regions. 

Using the differential partonic cross sections reported in the previous sections, one 
finds the explicit 2;*-dependencies of ^ce for the three cases: 



^CE — ^CE 
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- z' [z' 
2 



*2 



3) 



A' 



G 
CE 



^6E + e 



G 



2z 



*5 



+ e 



G 



^CE 



e^^E + <(2^*-l). 



(3.2) 

(3.3) 
(3.4) 



In Eq. (EJD, e^, and e^^^ are the fractions of resonant events from the processes 



G l^l and from the SM background, respectively, with e^^ + e 
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1. 
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Analogous definitions hold for Eq. (13.41) . where now + e^'^ = 1. Their dependence 
on the dilepton invariant mass M is determined by the overlap of parton distribution 
functions in Eqs. (12.1 p and (12.21) . Actually, the above equations strictly hold for ^^ut = 1, 
while all the results and figures reported here will be obtained by taking detector cuts into 
account. Differences turn out to be appreciable, and have an impact on the assessment 
of identification reaches, only near 2* = 1, whereas the 'optimal' values used in the Ace 
analysis will be z* ^ 0.5. Thus, Eqs. (I3.2I) - (I3.4I) are adequate for illustrative purposes, 
while giving results essentially equivalent to the 'full' calculation. 

As shown by Eq. (13.21) . the peculiar property of the observable Ace, as a function of z* , 
is that it is the same for all spin-1 exchanges, the SM 7, Z and any Z' model, regardless 
of the actual values of the left- and right-handed coupling constants to fermions, of the Z' 
mass Mz' and, to a large extent, of the choice of parton distribution functions. Deviations 
of Ace from the SM predictions can therefore be attributed to spin-2 or spin-0 exchanges 
in ([EID. 

To assess the domains in which the spin-0 and spin-2 hypotheses can be excluded 
as sources of a peak observed at M = Mr-, while giving the same number of events in 
the assigned interval AM in Eq. (12. 3p we start from the assumption that spin-1 is the 
'true' origin of the resonance. The level at which the two alternative hypotheses can be 
excluded for all 'allowed' values of their relevant model parameters, is determined by the 
experimental data, from the prediction of spin-1 iV) exchange, from the predicted spin-0 
{S) and spin-2 {G) exchanges, respectively: 

^^CE = ^ce ~ ^CE ^Aqe = — Aqe- (3.5) 

Of course, the identification potential will depend on the available statistics, i.e., from the 
number of signal events collected in the assigned M-interval, in addition to the systematic 
uncertainties. The latter are however expected to have a reduced influence on Ace, 
because it is a ratio of cross sections. 

As an example in Fig. O left panel, the center-edge asymmetry Ace is depicted as a 
function of z* for resonances with different spins, the same mass Mr = 3 TeV and the 
same number Ns of signal events under the peak. As anticipated, the calculations are 
performed using detector cuts and, also, the SM background has been accounted for. The 
deviations (13.50 are plotted in the right panel of the figure. The vertical bars attached 
to the solid line represent, again as an example, the 1-cr statistical uncertainty on the 
Aq-^ corresponding to the Z^-model with the assigned mass Mz'- Fig. [2] shows that the 
Z^-model with mass Mz' = 3 TeV can be discriminated from the other spin-hypotheses 
at the 2-a level by means of Ace at z* ~ 0.5. 

While A^g is independent of energy, Riid. ^q;^ are not. In the limit of little 
background, e^''^ will be small, and A^^ will only depend weakly on the energy. On 
the other hand, even in this limit, A^^ will in general have a significant dependence on 
energy, via the relative magnitudes of the fractions and e^. An exception to this energy 
dependence is the region around z* ~ 0.5, where the coefficient of vanishes. In this 
case, we have 

A^^iz* ^ 0.5) > A'c^iz* ^ 0.5). (3.6) 

This property is of course reproduced in Fig. [21 and allows to conclude that, in order to 
identify the spin-1 Z' resonance, if one is able to exclude the spin-0 hypothesis, the spin-2 
graviton of the RS model will then automatically be excluded. 
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Figure 2: Left panel: Aqe vs. z* for resonances u (spin-0), Z' (spin-1) and G (spin-2) 
with equal masses of 3 TeV. The error bars are the statistical uncertainties at the 1-a 
level on A^g"'^ for the ■i/'-model at 100 fb~^. Right panel: Asymmetry deviations, IAAceI, 
of the spin-0 and spin-2 hypotheses from the spin-1 one, compared with the uncertainties 

on /±r-p, . 



In order to determine the spin-1 signature space where the spin-0 hypothesis could 
be excluded against the spin-1 one, the deviation fl3.5p should be compared with the 
statistical uncertainty on Aqe expressed in terms of the desired number (k) of standard 
deviations {k^ = 3.84 for 95% C.L.). Notice that in practice Ace is almost unaffected by 
systematic uncertainty being a relative quantity. We have the condition 

\AA^^\ = k-6A^^, (3.7) 

where, taking into account that numerically (Aq^)^ ^ 1 at 2* ^ 0.5, 



*^c^E = y^^P^«^^. (3.8) 

V ' mm V ' mm 

From Eqs. fl3.7p and fl3.8p . one therefore obtains 

-^min — -^min ~ ( y aV ] ' (3-9) 

using the fact that — 0-5)1 is small compared to |A((g(2* ~ 0.5) |, and AA^^ = 

—s^Aq^ with = 1 — ^ 1. From Eq. (13. 9p one can easily evaluate the minimal 
number of events required to exclude the spin-0 hypothesis and, automatically, spin-2 as 
well, and in this way to establish the spin-1. One finds using Eq. (13.91) for the exclusion 
of the spin-0 resonance at 95% C.L., A'min — 110. One should emphasize that A^min 
determined from Eq. (13.91) is a model-independent value, since A^E defined in Eq. (ESD 
is independent of specific Z' models, being 'universal' for all spin-1 intermediate states. 
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Accordingly, spin-1 of the discovered resonance can be established at 95% CL. if resonance 
event samples Ns at the level of A^min or larger would be available. 

The behavior of N^^^ vs. Mr, as presented in Fig. [H is derived from the full calcula- 
tion including detector cuts, using the general Eq. (13. 7p . The intersection of the curves 
describing Ns against Mr for specific Z' models and displayed in Fig. [1] with the line of 
Amin determines the values of the Z' masses where the spin-1 hypothesis can be identified. 
One finds that for Mz' < 3 TeV the spin of Z' can be determined at 95% C.L. for all 
models under study, at 14 TeV, with Ci^t = 100 fb~^. For completeness we also display 
in Fig. [T] A'^i^ that lies below A^^^in anticipated. 

In addition to the illustrative consideration above, one can quantify the identification 
reach on the spin-1 hypothesis performing a 'conventional' analysis to obtain the 
exclusion domains of the spin-2 and spin-0 hypotheses. In this case, the function is 
defined as: 



AA 



CE 



A4 ' (3.10) 

O^CE 

with AAqe the deviations (13. 5p and 6Ace the statistical uncertainty (a specific spin-1 Z' 
model is taken the 'true' one) 

SAcE = a/ V^^H - (3-11) 




Figure 3: Discovery limits on Mz' (5-0" level) and Z'-spin identification reaches (95% 
C.L.) for neutral gauge bosons of representative models, using the lepton-pair production 
cross section a ■ Bi (/ = e,/i) and center-edge asymmetry Ace, respectively, at the LHC 
with integrated luminosity of 100 fb~^. Also, Z'-model distinction reaches (95% C.L.) are 
obtained from the analysis of the leptonic event rates. 

Like before, the spin-1 V model can be assumed to be the 'true' one, and the (95% 
C.L.) exclusion domains of spin-0 and spin-2 can be determined. Here, we combine the 
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channels I = e,fi. The 95% CL. identification reach of the spin-1 hypothesis in the 
signature space then results from the domain complementary to the combination of the 
spin-0 and spin-2 exclusion domains. The results of this numerical analysis are represented 
in the signature space (M^/, Ns) in Fig. [TJ In fact, the spin-0 exclusion is more restrictive 
than that for spin-2, as discussed above. Fig. [3] gives the 'translation' of the discovery 
reach on Z' models as well as identification reach on Z' spin presented in Fig. [T], in the 
form of a histogram. As one can see from Fig. [HI the spin-1 identification (or, actually, the 
spin-0 and spin-2 exclusion) can be obtained up to Z' mass of the order of 2.5-3.5 TeV, 
depending on the specific models. The model dependence of the spin identification reach 
is due to the difference in statistics, stemming from the different cross sections associated 
with these models. 

It might be useful to conclude this section by emphasizing that the basic observable 
Ace of Eq. fl3.ip only uses the z = cos ^c.m. even part of the angular differential distribu- 
tion, similar to the total cross section. The 'center' and 'edge' angular integration regions 
are symmetric around z = and include events with both signs of z. This might mitigate 
the impact of the ambiguity, at the pp colliders, in the experimental determination of the 
sign of z that can affect observables sensitive to the 2-odd part of the angular distribution. 

4 Differentiating Z' models 

Once the spin-1 character of a Z' peak at M = Mr = Mz' has been verified by the ex- 
clusion of the spin-0 and spin-2 hypotheses, one can attempt the more ambitious task of 
identifying the resonance with one of the Z' models by means of the measured production 
cross section a ■ Bi or, equivalently, of the peak event rate Ns- As anticipated in Sub- 
sec. 12.41 Bi will be assessed under the simplifying assumption of Z' decays to SM fermions 
only. Results from other observables, such as a ■ Tz', Aps, etc., have been qualitatively 
summarized in Sec. [T]0 

One can see from Fig. HI representing Z' signal event rates for two specific values of 
the Z' mass, that at Mz' = 2 TeV, an integrated luminosity of 100 fb~^ will be sufficient 
to distinguish all considered models pairwise, whereas at Mz> = 3 TeV, one is unable to 
distinguish neither between Z^ vs. nor between Z^ vs. Zlr. 

To perform this discrimination, we make the hypothesis that one of the considered 
models (ZgsM' ^i>y ^'rjy ^x' "^LR' ^alr) is the 'true' one, compatible with the measured 
cross section, and test this assumption against the remaining five models that in general 
can predict, for the same Mz', a different value of the cross section but within the uncer- 
tainty band of the former, hence not distinguishable from it. As a simple criterion, one 
can define a 'separation', in peak event rates Ns, between the 'true model' and the others, 
and then associate the foreseeable identification reach on the chosen 'true' model to the 
maximum value of Mz> for which all five such separations are larger than an amount 
specified by a desired confidence level. Finally, one can iterate this numerical procedure, 
in turn, for all six considered Z' models. 

For definiteness, we work out explicitly the example of the identification reach on the 
Z^LR model. We introduce the relative deviations of the event rates predicted by this 

^Actually, a precise measurement of the ratio Tz' /Mz' , if feasible, might also represent a discrimination 
criterion among classes of Z' models by itself. 
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Figure 4: Resonance event rates obtained for the reference models at Mz' = 2 TeV (left 
panel) and 3 TeV (right panel) are displayed. The error bars show the 1-a statistical 
uncertainty at 100 fb~^ of integrated luminosity. 



model, at the generic M = Mz', from those predicted by the other Z' models: 



ANs 
Ns 



Ns{Z') 



Ns{Z'alr) 



(4.1) 



Figure [5] shows the relative deviations (14. ip as functions oi Mz'- Vertical bars represent the 
1-a combination of the statistical uncertainty predicted by the ALR model at integrated 
LHC luminosity C^t = 100 fb~^, with the major systematic uncertainty for the total cross 
section a-Bi, represented by the uncertainties on the parton distribution functions (PDF). 
These are calculated using the CTEQ6.5 NLO PDF sets [21]. 

Corresponding to the definition (14. ip . one can introduce a function 



X 



ANs 
6N, 



(4.2) 



with SNs the corresponding experimental uncertainty, which includes both the statistical 
and systematic errors combined in quadrature, the former being determined by the .^alr 
model prediction of the event rate. It turns out that for M ^ 4 TeV the systematic 
uncertainty is larger than the statistical one, they cross over at a value around 15%. This 
systematic uncertainty has an effect comparable to that of the statistical uncertainty for 
Z' resonances with masses larger than 3 TeV [7]. 

On the basis of such ^6 can determine the maximum value of Mz' (hence the 
identification reach) for which the Z' models, with Z' ^ -^alr; t)e excluded once the 
ALR model has been assumed to be the 'true' one. This value must satisfy the conditions 

> XcL ^'^5 where Xc l determines the chosen confidence level. The results of 

this procedure, applied in turn to all six Z' models, are reported in Fig. [3|, together with 
the discovery and the spin-1 identification reaches, and indicate that all models under 
considerations might be distinguishable up to of the order of 2.1 TeV. 

Of course, these results rely numerically on the assumption about Z' decay stated 
at the beginning. On the other hand, neither of the Z' curves in Fig. [1] intersects with 



13 
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(TeV) 



Figure 5: Absolute value of relative deviation of the number of events, Eq. (14.11) . from 
the ^ALR predictions shown in Fig. 1 as a function of M^/ for neutral gauge bosons of 
representative models. The error bars are the uncertainties at the 1-a level on ANs/Ns 
for the ALR model, C.^t = 100bf~\ 

any other, so that the simple (and directly measurable) total cross section might also be 
considered a natural discriminator among models. 

5 The reduced-energy, low-luminosity case 

It may take quite some time before the experiments will be able to collect 100 fb~^ of data 
at 14 TeV. In the present section, we indicate how the spin-identification reach depends 
on the integrated luminosity, at two different energies, 10 and 14 TeV. The former is 
expected to be the energy initially available at the LHC 

Of course, in general, a detailed assessment of the spin-identification dependence vs. 
the total energy and the Z' mass depends on the, in some cases complicated behaviour of 
the individual parton distribution function, in addition to the applied cuts. Nevertheless, 
for a simplified estimate leading to an explicit, parametric, expression of such energy de- 
pendence, in a very rough sense one can trade energy for luminosity. Indeed, according to 
Sec. [3], the procedure of Z' spin-1 determination using Ace basically consists in excluding 
the spin-0 case, since spin-2 is then automatically excluded. The spin-0 cross section is 
determined by the parton cross section together with the overlap of the d and d parton 
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distribution functions in the proton, 



J,j(s, M) = J dym,, M)/,-(e2, M). (5.1) 



The spin-1 cross section also depends on the corresponding overlap of the u and u parton 
distribution functions, /^^(s, M). In the approximation that these two overlap integrals 
have the same dependence on E = ^/s, and in the narrow-width approximation, all these 
cross sections would have the same energy dependence, but they would differ by constant 
ratios. Thus, the spin-identification reach would be propoprtional to a unique function of 
energy, or equivalently, integrated luminosity, the same function for all Z' models. This 
reach would be given by the number of events, proportional to the integrated luminosity, 
the parton cross section a, and the overlap integral: 

N C,,,,aI{s,M). (5.2) 

The spin identification reach basically requires a certain number of events (around 110 
at the higher values of Mz'), i.e., the scaling with energy and luminosity is determined 
by keeping the above expression fixed as the luminosity or energy is changed. In the 
exponential approximation to the overlap integral, given by Eq. (3.16) of Leike [12] (valid 
for a wide range of E/Mz'), this becomes 



~ Ant^exp 



C2MZ' 
E 



(5.3) 



Neglecting the energy and mass dependences of Ci and C2, one finds a reach in Mz' that 
grows linearly with log£int, with a slope proportional to the beam energy, features which 
are qualitatively reflected in Fig. El which displays the spin-1 identification reach vs. the 
Z' mass for the various models at the two energies mentioned above. Experimental cuts 
have been taken the same as detailed in Sec. [2] for both cases. 

Figure [U] speaks for itself. It shows that at 10 TeV, and depending on the model, with 
a considerable fraction of 1 fb~ one could in principle identify the spin-1 of a hypothetical 
Z' with mass in the range 1-1.5 TeV. This range is expected to be in the reach of the 
Tevatron, for the planned luminosity increase from 2.5 up to 9 fb~^ [23]. However, the 
spin identification would be a unique feature of the LHC. 



6 Summary 

We can summarize the main part of this paper, relevant to the Z' identification at 14 TeV 
and 100 fb~^ luminosity, as follows: if new heavy resonance peaks will be discovered in 
the dilepton mass distributions for process (11.11) . a Z' can be observed up to Mz> ~ 4 — 5 
TeV. The statistical significance of measurements of the evenly-integrated (in cos^c.m.) 
asymmetry Aqe will allow to establish the spin-1 (or, to exclude the spin-0 and spin-2) 
of a heavy Z' gauge boson for Mz' < 3.0 — 3.8 TeV, at 95% C.L. We also assess the 
mass limits on Z' for which the studied Z' models can be distinguished, besides their 
common spin-1, in pairwise comparisons with each other. By a simple criterion based on 
the expected statistics, we find that with 100 fb~^ of integrated luminosity one should be 
able to distinguish among the six Z' models up to Mz' — 2.1 TeV (95% C.L.). 
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Figure 6: Spin-determination reach as a function of integrated luminosity, for 10 TeV 
(dashed) and 14 TeV (sohd). 
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